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Foreword 


a.  The  objective  of  this  program  was  to  develop  practical 
design  techniques  for  broadband,  phase  stable,  high  power, 
klystron  microwave  amplifiers  in  a  concise  mathematical  form. 

The  analytic  model  developed  was  to  be  accurate  at  all  signal 
levels  and  sufficiently  general  to  include  significant  space  charge 
modes.  The  predicted  interaction  was  to  be  investigated  for 
synthesis  methods  for  the  design  of  amplifiers  with  specific  charac¬ 
teristics.  Sample  design  calculations  for  an  aunplifier  with  supplied 
parameters  was  to  be  carried  out  following  the  successful  develop¬ 
ment  of  the  theoretical  model. 

b.  The  contractor  developed  a  statement  of  the  beam  space 
charge  interaction  in  the  linear  or  small  signal  section  of  the 
amplifier  in  matrix  notation  and  from  this  deduced  a  gain  equation 
in  closed  form  for  the  synchronously  tuned  cavity  case.  The  re¬ 
sulting  amplifier  design  required  many  cavities  tuned  to  a  frequency 
above  the  operating  frequency.  This,  however,  produced  a  very 
high  gain  peak  at  the  cavity  resonant  frequency  which  required  elim¬ 
ination  to  prevent  oscillation.  The  calculations  necessary  to  slightly 
detune  the  cavities  and  eliminate  the  gain  peak  were  not  reducible 

to  any  short  form  so  a  computer  would  be  required  to  make  the  cal¬ 
culations.  The  large  signal  beam  interaction  theory  was  developed 
in  the  infinite  beam  model  for  one  gap  and  drift  tube,  however,  the 
initial  conditions  encountered  at  the  second  gap  precluded  a  solution 
in  closed  form  for  the  two  or  more  gap  case.  The  sample  design 
was  completed  in  a  very  elementary  form  because  of  the  extremely 
short  time  remaining  in  the  contract  period  when  the  theoretical 
portion  of  the  program  was  completed. 

c.  The  design  techniques  developed  in  this  contract  have  been 
developed  in  theory  only.  Experimental  klystrons  developed  by 
other  contractors  have  demonstrated  the  validity  of  some  of  the 
results  predicted  by  this  theory  but  no  amplifier  has  been  built  using 
this  design  technique.  This  large  signal  theory  has  eliminated 
several  important  approximations  made  in  other  theories  but  requires 
more  work  to  be  completely  developed  for  two  or  more  gap  interactions 

d.  An  analytic  model  of  the  gain  equation  for  a  klystron  type 
amplifier  which  allows  synthesis  of  amplifiers  with  broadband  and 
controlled  frequency  characteristics  will  be  of  considerable  use  in 
the  design  of  amplifiers  for  advanced  radar  systems. 


Frank  E.  Welker,  1/Lt,  USAF 

Project  Engineer 

RADC,  Griffis s  AFB,  NY 
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Abstract 


Matrix  methods  are  used  in  an  attempt  to  clarify  the  predic  - 
tions  of  many  resonator  klystron  response  at  frequencies 
above  the  pass  band.  The  gain  expressions  have  been  ex¬ 
pressed  in  factored  form  in  the  case  of  staggered  tuning  but 
numerical  calculations  in  the  super  pass  band  region  have 
not  been  completed  due  to  computational  length  and  complexity. 

Several  subjects  in  large -signal  beam  theory  are  treated.  The 
multigap  problem  is  formulated  in  terms  of  an  integral  equation 
which  is  solved  within  the  framework  of  second  order  pertur¬ 
bation  theory.  The  disc -model  of  space -charge  forces  is  used 
to  derive  a  large -signal  equation  of  motion  for  a  finite  diameter 
beam.  Solutions  to  the  equation  indicate  different  space  charge 
reduction  factors  for  each  harmonic  of  the  modulating  voltage. 

Statistical  mechanical  methods  are  used  to  derive  equations  of 
motion  corresponding  to  the  infinite  beam  model.  Terms  ap¬ 
pear  which  cast  doubt  on  the  infinite  beam  model  near  cross¬ 
over.  The  new  terms  are  simulated  by  the  introduction  of 
phenomenological  collision  loss  terms  which  impede  the  onset 
of  crossover. 

Our  results  still  generally  uphold  the  idea  that  the  klystron 
can  be  built  as  a  broadband  device  with  sophisticated  phase 
versus  frequency  characteristics. 
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I.  General  Introduction 


In  an  earlier  engineering  report^  there  appeared  among 
other  things  an  analysis  of  the  gain  eind  phase  versus  frequency 
characteristic  for  a  many -resonator  klystron.  The  analysis 
was  based  upon  a  single  space  charge  mode -small  signal  model. 

In  addition,  empirical  evidence  was  cited  as  a  basis  for  pre¬ 
diction  of  performance  of  the  high  signal  level  gaps  and  drift 
tubes  for  a  very  broad  band  klystron  having  excellent  phase 
linearity  in-band  and  normal  klystron  efficiency  throughout  that 
band. 

The  analysis  consisted  of  manipulations  of  well-known 
space  charge  relationships.  The  resulting  design  called  for 
tuning  of  the  many  resonators  to  resonant  frequencies  outside 
of  and  above  the  pass  band  frequencies.  This  produced  a  gain 
characteristic  flat  and  linear  in  phase  within  the  pass  band  but 
it  resulted  in  the  possibility  of  an  enormously  large  gain  at  one 
frequency  above  the  pass  band.  Stagger  tuning  was  recom¬ 
mended  as  a  means  of  eliminating  the  large  responses  above 
the  pass  band.  While  the  analysis  is  believed  to  be  most  con¬ 
vincing  regarding  the  pass  band  frequencies,  it  leaves  some¬ 
thing  to  be  desired  in  terms  of  conviction  at  the  higher 
frequencies. 

1.  RADC-TDR-62-199,  J.  F.  Kane  and  R.  N.  Wilson, 
"Advanced  Klystron  Study,"  Kane  Engineering  Laboratories,  Palo 
Alto,  Calif.,  prepared  under  Contract  AF30(602) -2423  (June  1962) 
AD-278  223. 
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Much  of  the  work  being  reported  here  was  aimed  at  clari¬ 
fying  the  predictions  of  tube  response  at  frequencies  above  the 
band.  Exact  expressions  have  been  obtained  for  this  region  as 
well  as  the  pass  band  region  but  calculations  remain  laborious  in 
this  region  despite  very  considerable  efforts  to  find  means  for 
simplifying  them.  However,  calculation  of  pass  band  response  is 
a  relatively  simple  matter  for  those  cases  which  call  for  great 
phase  linearity  in  the  pass  band. 

We  have  simplified  the  derivations  of  formulae  here.  We 
find  it  necessary  in  this  report  to  resort  to  a  detailed  calculation 
of  the  tube's  response  at  the  frequencies  above  the  pass  band  as  a 
means  of  revealing  the  efficacy  of  stagger  tuning  for  annihilation  of 
the  super  pass  band  response  peak.  The  numerical  calculation  of 
gain  in  the  super  pass  band  region  has  not  been  completed  due  to 
computational  length  and  complexity. 

Our  work  in  large -signal  theory  has  progressed  in  several 
areas.  We  have  formulated  the  multigap  problem  using  the  infinite 
beam  model  for  space -charge  forces.  We  obtain  an  integral  equa¬ 
tion  for  the  displacement  vector  which  provides  a  convenient  starting 
point  for  numerical  calculations.  We  have  not  been  able  to  obtain 
closed  form  large -signal  solutions  to  the  multigap  problem  because 
of  the  difficulty  with  boundary  conditions  at  the  second,  third,  etc.  , 
gaps.  However,  we  have  solved  the  integral  equation  for  the  case  of 
two  gaps,  within  the  framework  of  second  order  perturbation  theory 
and  this  solution  displays  several  interesting  features. 


We  have  formulated  the  large -signal  problem  using  a 
disc -model  for  the  space -charge  forces.  The  equation  of 
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motion  takes  a  more  complicated  form  than  that  of  the  infinite 
beam  model.  We  content  ourselves  with  an  approximate  large - 
signal  equation  of  motion  which  we  can  solve  exactly  for  certain 
forms  of  gap  modulating  signals.  The  solutions  exhibit  space - 
charge  reduction  properties  and  it  is  shown  that  a  different 
space -charge  reduction  factor  occurs  for  each  harmonic  of  the 
modulating  voltage. 

We  next  develop  the  infinite  beam  model  starting  from 
statistical  mechanical  foundations.  The  resulting  equations  of 
motion  contain  heat  diffusion  terms  in  addition  to  the  usual 
terms  which  occur  in  the  hydrodynamic  equations  of  beam  theory. 
The  extra  terms  appear  to  become  large  at  crossover  and  this 
indicates  that  the  infinite  beam  model,  which  pictures  the  beam 
as  a  set  of  charge  planes,  may  cease  to  be  valid  near  crossover. 

We  take  a  step  toward  representing  the  diffusion  terms, 
mentioned  above,  by  introducing  a  phenomenological  collision 
loss  term  into  the  infinite -beam  model  equation  of  motion.  We 
are  able  to  solve  the  resulting  equation  of  motion  in  closed  form 
in  the  case  of  modulation  at  an  infinitesimally  extended  buncher 
gap.  We  find  that  the  presence  of  the  loss  tends  to  impede 
crossover. 

Finally,  we  remark  on  some  of  the  directions  which  future 


large -signal  beam  theory  might  take.  There  are  still  many 
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interesting  and  difficult  problems  in  large -signal  beam  theory  which 
remain  to  be  solved.  In  fact  we  may  claim,  so  far,  only  to  have 
scratched  the  surface  in  this  particular  branch  of  the  broad  field  of 


non-linear  oscillations. 


-  5  - 


II.  N  Resonator  Gain  Formula 

We  identify  a  drift  tube  by  assigning  to  it  the  number 
of  the  resonator  which  precedes  it.  The  resonators  are 
numbered  consecutively  from  1  to  N  starting  at  resonator 
number  1,  the  input  resonator.  We  are  interested  in  the 
peak  values  of  the  a-c  current  and  velocity  modulation  on 
the  electron  beam.  These  are  denoted  as  follows. 

a -c  current  entering  the  n^^  resonator  i^' 

til 

a-c  current  leaving  the  n*""  resonator  i^ 

a-c  velocity  entering  the  n*^  resonator  v^' 
a-c  velocity  leaving  the  n^^  resonator  v^ 

In  order  to  work  with  dimensionless  matrix  elements 
we  shall  measure  a-c  velocities  in  current  units  as  given 
below. 

For  the  n^^  gap  the  following  transformation  applies. 
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where  in  MKS  units 


g 


z 

n 


M  —  Z  (oj)  ,  gz  IS  dimensionless  J 
n  u  n  ®  n 

o 

is  an  angular  measure  of  time  delay  in  the  gap. 


(2) 


CO 

a 

e 

o 

I 

o 

u 

o 

r 


M 

n 


Z  (w) 
n  ' 


is  the  operating  angular  frequency 

is  the  electron  beam's  radius 

is  the  dielectric  coefficient  of  empty  space 

is  the  total  d-c  beam  current 

is  the  specific  electronic  charge 
is  the  d-c  beam  velocity 

is  the  plasma  reduction  factor  for  the  beam 
and  drift  tube 

is  the  gap  coupling  coefficient 

is  the  gap  impedance  at  the  operating  frequency 


4^  measures  transit  time  delay  in  radians  of  the 
operating  frequency  in  the  structure  to  which  the  transfor¬ 
mation  applies.  Due  to  the  linearity  of  the  formalism,  only 
the  total  phase  delay  (flight  time  for  beam  electrons  traversing 
the  entire  tube  structure)  is  important.  Therefore  we  have  not 
bothered  to  attach  a  phase  delay  to  each  tube  section.  Rather, 
the  total  delay  is  prefixed  to  the  composite  tube  gain  expression 
at  a  later  point  in  the  analysis. 
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For  the 


drift  tube  the  following  transformation  applies. 


j  \  / 4h- 1 


(3) 


where  6  =  {cj  d  )/u  :  u  =  ru  euid  u  is  the  beam  plasma 

n'qn  o  q  p  p  ^ 

angular  frequency,  d^  is  the  length  of  the  drift  tube. 

We  a.pply  the  drift  tube  and  gap  transformations  in  the 
same  order  as  these  are  encountered  by  the  electron  beam  in 
traversing  the  tube.  The  first  resonator  gap  calls  for  special 
attention  because  it  is  externally  driven  and  is  supplied  with 
an  unmodulated  becun  in  all  cases  where  the  drive  power 
supplied  is  large  compared  to  available  noise  energy  on  the 
entering  beemi. 


To  account  for  first  gap  conditions  we  set 


This  is  equivalent  to  the  assumption  that  a  peaJc  r-f  voltage 

exists  on  the  first  gap  at  all  frequencies.  Thus  the  circuits 
used  in  applying  external  drive  have  been  temporarily  excluded 
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as  factors  in  the  broadband  properties  of  the  rest  of  the  tube. 
Now  entering  the  n^^  resonator  we  have 


where  i|i.  -  wd,  /u  and  d,  is  the  distance  between  the 
t  In  o  In 

th 

centers  of  the  1st  and  n  gaps. 

th 

Noting  that  the  peak  r-f  voltage  on  the  n  gap  is  given 
by 


I  / 


Z.  CJ) 


(6) 


we  can  write 


(7) 
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where  K  = 


gM,-!l-M  Z  (w) 
*  lu  n  n'  ' 


V^/Vj  is  of  course  the  gain. 


o 

The  relatively  complicated  quantity  representing  the  modified 
velocity  is  of  no  immediate  interest  in  a  gain  calculation. 


The  matrix  (D^)  has  eigenvalues 
diagonalized  by  U  where 


and  is 


Vrx  V.  \ 

'In  -‘/rj 


(8) 


with 


l)=  0"'  ,  (JU''=  I 

A 

We  write  (G)  in  the  form 


(9) 


and  make  use  of  the  following  definitions. 
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UtPnW 

\0 


\ 

Lents, J 


U(.G-k)0  = 


(12) 


14  Vh 

-  V*,  H 


(13) 


(14) 


-  ll  - 


A 

The  vector  (G)  now  becomes 


(15) 


If  we  apecialice  the  tube  design  by  setting  all  n  drift  tube 


lengths  equal  we  have 


(16) 


We  notice  that 

(©n)^0)"  (Sj  +  (l) 

where 

Cl)  =  ("!) 

Co^  M  =  9  +  J'l.  j  5  ^  ^ 


(17) 


(18) 


This  is  the  same  notation  that  was  used  in  the  previous 

♦  (2) 

engineering  report. 

In  addition,  we  observe  that 

©P^  U  j  *  (^) 


2.  RADC-TDR-62-199. 
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Multiplying  the  first  few  matrices  on  the  right  of  the  gain  ex 
pression  into  the  vector  6  we  get 


©  (\  ©  =  Ux  {s)  +-  Cl ) 


(20) 


•— 1  w  ^  © 

I 

+  (2c<>^M3  1- 


(21) 


(22) 


etc.  Observation  of  the  above  polynomials  and  the  manner 
in  which  they  are  generated  shows  that  the  gain  can  be  written 


~  ^  K  0  I  ^N-i  1 


(23) 


where 


•XCiK-Ctti.  { 


I 


o 


(24) 


o 


I 


which  determinant  is  (N  -  2)  by  (N  -  2). 

If  there  is  no  stagger  tuning  the  are  identical  and 
the  gain  becomes 


This  result  we  have  seen  before. 
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III.  Computation 

We  have  searched  for  simple  means  for  evaluating  the 
determinant  )^u|  such  that  a  tube  design  engineer  could  per¬ 
form  design  calculations  at  a  desk  in  a  few  hours.  We  have 
not  as  yet  succeeded  in  this.  To  perform  the  calculations  for 
a  tube  of  fcur  or  more  low  level  resonators  is  still  a  very  labo 
rious  undertaking  without  a  computer.  Computer  programs 
exist,  however,  which  would  permit  stagger  tuned  pass  band 
computations  in  a  few  minutes.  In  consequence  of  this,  we 
have  developed  the  following  technique  for  desk  computation. 

The  gain  determinant  ,  when  written  out,  appears 

as  the  quotient  of  two  polynomials  in  the  modulating  frequency 
We  wish  to  express  these  two  polynomials  in  factored  form 
for  the  case  of  non -synchronous  tuning.  Since  the  factor¬ 
ization  is  somewhat  involved  we  first  present  the  general 
results.  It  is  to  be  noted  that  the  factorization  is  performed 
separately  for  the  numerator  and  denominator,  thus  avoiding 
the  pitfalls  of  perturbation  theory  in  its  usual  sense 

At  the  end  we  present  a  few  special  cases. 
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Results 


We  express  the  gain  determinant,  given  by  Eq.  (39) 
below,  in  factored  form  when  the  cavities  have  resonant 
frequencies 


«  I 


(26) 


and  we  allow  a  change  in  drift  tube  lengths 


€>  - ^  €>,  +  A  ^  , 


(27) 


For  an  (N  -  1) -dimensional  determinant  we  obtain 


where 


(28) 


= 


+  1+iId,+ 


C.M- 


(29) 


(30) 
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N  =  K  a  (cj  -  a  )  (w  -  b  ) 
o  o  o  o  o 


N.  =  -K^q^(u  -  a^)  (w  -  b^)  (c  +  tan  0  A  0) 


J 


D  =  (bj  -  c  )  (oj  -  d  ) 
o'  o'  o' 


D.  =  cj  (d  -  c  )6. 
j  o  o  o  j 


K  =2  cos  0 
o  o 


tan  0 


b  =  j 


.G 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 

(37) 


(*>oi  .  Qq>  Qq  are  given  by  Eq.  (48)  below. 

Factorization 

We  begin  by  rewriting  the  gain  determinant,  given 
by  Eq.  (24),  in  the  following  form.  Note  that  N  of  this 
section  is  to  be  identified  with  N  -  1  of  Section  II. 
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(39) 


This  is  accomplished  by  some  interchanging  of  rows  and 
columns  in  Eq.  (24).  We  examine  cosh  u 


(40) 


We  wish  to  express  A  in  factored  form  when  Uj^  4 
4  4  •  •  •  4  j  .  We  assume  that  the  only  frequency 

dependence  of  g  is  in  the  multiplicative  w  term.  Define 


(41) 


where  D  is  a  constant  frequency  chosen  equal  to  the 
resonant  frequency  of  the  cavities  in  the  synchronous  case. 


(aJ 


o 


\ 

Vc 


(42) 
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n  may  be  given  any  value.  It  occurs  only  for  dimensional 
reasons. 

Write  2  cosh  u^  in  factored  form .  The  synchronous, 
2  cosh  u  ,  is 


r.-i 

[i-iUal*-#] 

-Mix/  -U*  J 

1 

(43) 


as  can  easily  be  verified  from  Eq.  (40).  Introduce  the  notation 


\  —  t  to  —  J 


where 


K  =  2  cos  6 


1  -  tan  0 
2i2c 


J9 

c 


Factoring  Eq.  (43)  we  can  write 


(u  -c.)  (vo-  a.; 


(44) 


(45) 


(46) 


where 


(47) 
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f 

U  c 


UJ« 


(48) 


We  consider  the  case  where 


(49) 


Then 


(50) 


where  we  have  allowed  l/(4Qj  )  to  be  comparable  with  unity. 

To  treat  the  non -synchronous  case  we  start  with  the 
synchronous  case  and  consider  the  effect  of  changing  the 
resonant  frequency  of  the  cell.  We  also  allow  a  change 
in  drift  tube  length.  Eq.  (39)  is  only  valid  when  all  drift 
lengths  are  equal.  Hence  we  change  all  the  drift  tube  lengths 
by  the  same  amount.  We  let 


- >  U3*  Au>, 


nS'  0  C.  I } 


(51) 
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<<  » 


(52) 


0  - 00  +  ^0  ^  69  «\ 


(53) 


Then  keeping  only  lowest  order  terms 


(54) 


where 


2  AC- 


6o 


6  6 


I  - 


_iL 

lA  C 


"t 


I 


(56) 


Eqs.  (54),  (58)  tell  us  how  the  zeros  and  poles  of  Eq.  (46) 
move  when  we  change  the  resonant  frequency  of  the  cavity 
and  the  drift  tube  length  associated  with  cell  j  . 

To  get  the  perturbed  Z  cosh  u^  we  only  have  to  put 
in  the  changes  in  K,  a  . 

—  ko  ^  1“  ^  6 
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(57) 


(58) 


(59) 


(60) 


(61) 


oi.  ti-e) 
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With  Eqs.  (54),  (58),  (59),  (61)  we  obtain  2  cosh  u,  , 


(u-a«  -^ij) 


(62) 


Note  that  we  have  treated  the  movement  of  the  zeros  and  poles 
separately  and  they  are  constrained  in  their  movement  by 
Eqs.  (51),  (52),  (55),  (56),  (6 1)  in  order  that  the  laws  of 
physics  remain  satisfied.  That  is,  we  allow  only  move¬ 
ments  of  the  poles  and  zeros  which  are  physically  realizable. 
Whether  the  chosen  movements  are  practically  realizable  is 
another  question. 

Now  we  turn  to  the  problem  of  evaluating  the  gain  deter¬ 
minant  given  by  Eq.  (39).  First  we  note  the  following  properties 
of  a  special  matrix  and  its  determinant.  Consider  the  matrix 


The  eigenvalues  are 


A;  =  Y  + 

The  matrix  which  diagonalizes*"^^  is 


v>=  '  ^If  4,»>' 


We  now  write  Eq.  (39)  with  the  aid  of  Eq,  (62) 

*  US- 


where 


I*-'  . 

IT  lU -C,  -  J* +‘-*0^^) 


H,tW.  0,fP. 

p,+  0i  PcfOx 

O  +  Dj^^ 


A>\  = 


^  P»+D„_, 
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where 

D  +  D.  s  (u  -  c  -  u  6.)  (b)  -  d  +  w  6.)  (69) 

o  j  '  o  o  '  o  o  j'  '  ' 

N  -I-  N.  =  K  a  (1  -  «)  (  1  -  tane  A  e)  • 
o  j  o  o  '  '  ' 

.  (70) 

•  (<-  -  ao  -  A  ij)  («  -  bo  -  A  2j) 

Note  that  we  have  calculated  the  effect  of  tuning  separately 
for  the  numerator  and  denominator  of  .  That  is,  we 
calculate  the  movement  of  the  poles  and  seros  Separately 
and  in  this  sense  the  treatment  differs  from  perturbation 
theory.  Perturbation  theory  fails  when  we  must  manipulate 
resonant  quantities;  e.g. ,  move  poles  and  zeros  of  a  gain 
determinant. 

The  denominator  of  A  given  by  Eq.  (67)  appears  in 
factored  form.  It  remains  to  factor  the  numerator, 
given  by  Eq.  (68).  We  factor  A>|  through  first  order  in  the 
tuning.  Rewrite  77  in  the  form 


where 


(72) 
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Kq  (u-a)(u-b) 
o  o  o'  '  o' 


-K  a  (u  -  a  )  (u  -  b  )  (e  +  tan  6^6) 
o  o  o  o 

-K  a  [  A  , .  (w  -  b  )  +  A  , .  -  a  )] 

o  o  Ij  o'  2j  '  o'-* 


(u)  -  c  )  (u  -  d  ) 
'  o  o' 


(d  -  c  )6. 
o  o  o  j 


We  transform  0'\  with  the  matrix  U  given  by  Eq.  (65) 


and  note  that 


^  AO-*’) 


The  result  is 


a'h  = 


A  {yi',  +  ^-'^,  0} 

J 


where 
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X 


(77) 


N-l 

bX.-  A-K.'  -  T  {N.  +  lk ] 

^  =  1 


(V«) 


Also 


-  ....  ^ 

‘Hi.  ,  -  - - - 'Hlw-. 


V 


-  —  - - Aft-tVi*- 


(79) 


We  wish  to  calculate  through  first  order  in  the  quan¬ 

tities  .  To  see  what  this  is,  consider  expanding  tX  by 
minors  along  the  first  row. 

. 


(A.-+0 


"Hjv  ,Aj+i|,' ^ 


'■An., 


+  (Vitiii 
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Hjij  ----- 

***•  -----  '*l3[|»-» 

% 

I 

I 

%  ^ 

- -  Ji|l-l'*‘%-*,H-* 


(80) 


All  terms  in  (80)  except  the  first  are  second  order  in  the 
Ti'.j  and  therefore  are  negligible  in  our  approximation. 
Through  first  order  the  first  term  is 


hy\'r:  A 


I 


ij+Vj - 

♦  I 

'  I 

•  { 

•  I 

» 

- - •“ 


(81) 


Now  the  determinant  multiplying  in  Eq.  (81)  is  the  same 
form  as  our  original  determinant  AV  except  of  dimension 
(N  -  2)  X  (N  -  2).  Hence  through  first  order 
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'Hh . 


^11 


- - 'Im.w-* 


■  An -I 


4 


(82) 


"H jn  •- 


- ^3.^*/ 

► - '  ~'^1(,jvf 

i 

\ 

I 

I 

I 

I 

- 


(83) 


We  obtain  immediately  by  induction  that,  through  firit 
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order  in  the  tf* 


ij’ 


W-l 

=  II  A*  +  2^ 

j»'.  i*ix 


Eq.  (84)  y/e  express  in  the  form 

ATS  =  ■¥  . 


Thus  what  is  left  to  do  is  calculate 


(84) 


(85) 


(86) 


We  use  Eq.  (65)  and  (72)  to  obtain 

=  O.ilNiO.,  +D,  i>u) 


Hi>>a  +  PjW) 


+  0|,.,  (  D*-. 


(87) 
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Inserting  the  values  of  U^j ,  Eq.  (87)  becomes 

i*' 


Combining  Eq.  (85)  with  (88)  the  result  is 


From  Eqa.  (73),  (55)  and  (50) 
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"i” 


-  C***-**©)  t***"  ^»o)  (€•  A®  ) 


(91) 


-  ‘*3©  (  ^-o)  ^ j 


Substituting  in  Eq.  (90)  we  obtain 


H-l 

A'h  »  IT 


(92) 


where 

4't*^®.  A©)j(U>-a*);u»-lo) 


(u -Co)  ( 
N 


Eq.  (93)  is  a  quadratic  in  u  which  is  easily  factored. 


-  32  - 


Special  Cases 

We  look  at  the  following  special  cases. 

(1)  ^  • 

JAc 

We  consider  c  negligible  and  Eq.  (93)  becomes 

+  (U'k.) 

4  1  ^  C  ^  J (I- 

h 


In  this  approximation  the  gain  determinant  (39)  is 
independent  of  frequency  when  6j^  =  0  . 

Js-c,  ?= 


(95) 


-  ^k.-c.  (i  -  i^e.  AO)  + 


(96) 
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(2)  In  addition  to  the  conditions  of  (1)  we  let  ^  =  0  . 

Then 


h 


(97) 


(3)  Same  conditions  as  (2)  but  with  A  0  =  0  . 


1^0  T.  .lit 


4. 

+  li  i 


C  "I 

N 


The  special  case  6j^  =  (-1)  6  requires  N  -  1  be  an  even 
integer,  if 


=  o 


(98) 


(99) 


is  to  be  satisfied.  However,  we  find 
W-l 


kti 


in  this  case.  Hence,  to  the  order  of  our  approximation, 
this  type  of  stagger  tuning  has  no  effect  on  the  numerator 


(100) 
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of  the  gain  determinant.  It  makes  the  poles  of  order  (N  -  l)/2 
instead  of  (N  -  1). 


(4)  In  addition  to  the  conditions  of  (1)  we  let 


(101) 


K 

Hence 


(103) 


If,  in  addition,  A  9  =  0 


(104) 
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In  this  case 
M-l 


A  * 


W-  a,  (u  -k 


H-» 


(105) 


Hence  this  type  of  tuning  does  not  remove  the  confluence 
of  the  poles. 

From  (3)  and  (4)  it  is  clear  that  we  should  make 
all  6's  different.  For  example: 

(5)  In  addition  to  the  conditions  of  (1)  we  let  A  0  =  0 
and 


(106) 


This  will  make  all  poles  and  zeros  move  in  first  order. 
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IV.  Nvimerical  Computationa 

At  this  point  we  have  not  completed  any  numerical  com¬ 
putations  which  are  based  on  stagger  tuned  tubes.  As  a  result 
we  are  not  yet  in  a  position  to  demonstrate  by  actual  computation 
that  the  so-called  super  pass  band  gain  peak  can  be  annihiliated 
or  reduced  to  a  safely  small  maximum  value.  On  an  intuitive 
basis  we  do  not  question  that  this  gain  peak  can  be  annihiliated 
in  practical  structures.  For  the  design  of  broadband  structures 
with  gain  and  phase  linearity  as  required,  we  use  the  degen¬ 
erately  tuned  gain  formula  which  applies  in  the  region  below 
the  resonant  frequencies  of  the  resonators. 

Figure  1  shows  the  real  and  imaginary  parts  of  cosh  u  as 
a  function  of  frequency.  This  represents  a  tube  having  10  mega¬ 
watts  of  output  capability  using  a  beam  microperveance  of  unity 
near  2800  megacycles  per  second  for  which  6  =  0.  The  range 
of  frequencies  in  which  the  real  part  of  cosh  u  is  greater  than 
unity  and  within  which  the  imaginary  part  of  cosh  u  is  negligible 
is  of  the  greatest  interest  here.  So  long  as  the  imaginary  part 
of  cosh  u  is  small,  there  will  be  at  most  a  small  and  monotonic 
departure  of  the  phase  versus  frequency  curve  from  a  strictly 
linear  curve. 

Figure  2  shows  the  gains  whidh  can  be  expected  from 
klystrons  whose  small  signal  level  resonators  are  numerous 
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Fig.  1:  Real  and  imaginary  parte  of  coah  u 
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Fig.  2:  Gain  versus  frequency 

for  multiresonator  klystrons  in 
flat  phase  vs  frequency  region 
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and  closely  spaced.  The  drift  length  used  in  Fig.  2  is  2  inches 
so  it  can  be  seen  that  very  useful  tubes  can  be  built  in  lengths 
of  two  feet.  We  can  compare  this  with  tubes  which  provide 
45  db  of  gain  in  5  feet  of  length.  This  shows  the  dramatic  im¬ 
provement  available  in  terms  of  gain  per  unit  length  by  the  use 
of  numerous  closely  spaced  resonators.  It  can  also  be  seen 
in  Fig  1  that  there  is  a  near  zero  of  gain  at  a  frequency  of 
about  1-1/2%  above  the  resonant  frequency  of  each  cavity.  This 
null  would  be  so  located  in  stagger  tuning  as  to  radically  reduce 
in  size  the  gain  peak  which  would  be  obtained  in  the  case  of  de¬ 
generate  tuning. 

It  is  important  to  notice  how  slowly  the  real  part  of  cosh  u 
(in  Fig.  1)  sinks  down  through  the  value  unity  as  detuning  toward 
lower  frequencies  tadces  place.  There  is  less  than  10%  change 
in  going  from  5  =  .  1  to  -.2,  for  example.  Now  the  magnitude 
of  cosh  u  in  this  region  is  sensitive  to  the  value  of  perveance 
used  in  the  beam.  We  see  that  we  should  have  used  a  beam 
microperveance  of  about  2.  If  so,  we  would  have  obtained  a 
curve  showing  the  real  part  of  cosh  u  remaining  appreciably 
above  unity  for  more  than  15%  in  frequency  range.  Thus  we 
would  have  obtained  gain  magnitude  curves  such  as  those  of 
Fig.  2  which  are  higher  and  much  flatter  in  value.  In  addition 
we  would  have  found  ourselves  with  much  more  room  for  stagger 
tuning  of  the  resonators.  It  is  regretted  that  the  pressure  of 
time  has  prevented  our  carrying  out  such  computation. 
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In  a  tube  such  as  we  are  describing  it  is  important  that  the 
many  short  drift  tubes  present  no  feedback  possibilities  due  to 
drift  tube  propagation  at  harmonic  or  other  frequencies.  This 
is  especially  important  in  the  vicinity  of  the  resonant  frequencies 
where  gain  may  be  high.  Figure  3  shows  how  drift  tube  dimensions 


are  chosen. 
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Fig.  3:  Drift  tube  cut-off  wavelengths  and 

harmonics  of  the  operating  frequencies 


V.  High  Level  and  Output  Stages 


No  trustworthy  or  tractable  design  theory  exists  in  a 

form  useful  for  designing  the  high  level  klystron  stages.  In 

(31 

our  previous  report'  '  we  treated  these  by  reference  to  em¬ 
pirical  evidence. 

Also  in  our  last  report  we  developed  expressions  for 
a  multi -element  output  section  capable  of  presenting  to  the 
klystron  beam  an  optimum  value  of  output  gap  impedance 
over  an  adequate  band  of  frequencies.  We  find  no  need  to 
expand  upon  what  has  been  presented.  Some  concern  has 
been  expressed  to  the  author  regarding  the  bandwidth  pre¬ 
dicted.  This  concern  derived  from  comparison  of  predictions 
of  bandwidth  capability  with  similar  predictions  for  duplexer 
devices  which,  of  course,  must  be  matched  to  waveguide  im¬ 
pedances  throughout  the  band.  This  is  not  the  case  in  the 
output  of  a  modern  high-power  klystron,  however.  We  always 
desire  a  large  VSWR  seen  looking  back  into  such  a  klystron. 
Without  this,  tube  efficiency  would  either  be  severely  reduced 
in  one  case,  or  else  the  efficiency  would  be  cut  in  half  and 
excessive  r-f  voltage  would  appear  in  the  output  resonator 
in  the  other  case.  In  this  latter  case,  electrons  would  be 


3.  RADC-TDR-62-199. 
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reflected  back  through  the  drift  tube  in  a  manner  which  would 
tend  to  cause  spurious  oscillations  or  other  erratic  tube 
behavior. 

Some  brief  experimental  work  was  done  during  this 
program  toward  obtaining  measured  curves  of  broadband 
output  resonator  performance.  This  work  has  not  culminated 
in  presentable  data.  This  is  due  to  a  poor  quality  of  r-f 
contacts  used  in  the  cold  test  assembly.  It  was  considered 
to  be  not  of  sufficient  interest  to  undertake  the  expense  of 
preparing  soldered  output  sections  in  order  to  demonstrate 
what  is  believed  to  be  most  soundly  predicated  on  well-known 
information. 
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VI.  Infinite  Beam  Multigap  Theory 

1 .  Introduction  and  Summary 

We  present  here  the  results  of  some  work  on  the  modulation 
of  a  beam  by  many  gaps.  The  goal  was  to  obtain  a  large  signal 
description  of  the  multigap  problem.  The  results  obtained  have 
been  more  modest  due  to  difficulties  in  applying  boundary  con¬ 
ditions  at  the  2nd,  3rd,  etc. ,  gaps. 

The  infinite  beam  model  is  used  to  formulate  the  multigap 

(4) 

problem.  The  displacement  vector  formalism,  which  has  been 

( 5) 

successful  in  the  large  signal  description'  '  of  a  beam  modulated 
at  an  infinitesimally  extended  gap,  provides  us  with  an  integral 
equation  for  the  displacement  vector  in  the  multigap  situation. 

The  integral  equation  for  N-gaps  is  as  follows. 


(107) 


4.  RADC-TDR-62-199.  IV  Eulerian  Theory,  Section  A, 

p.  6l. 

5.  Ibid,  Section  C,  p.  91. 
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In  Eq.  (107),  o*  =  p  ^  is  the  dimensionless  displacement 
vector,  X  =  p  z ,  0  =  w/v  ,  a.  =  V.  /V  ,  6  =  ut ,  and  i. 
is  a  phase  angle  associated  with  the  i^^  gap.  This  equation 
provides  a  convenient  starting  point  for  numerical  calcu¬ 
lations.  We  have  used  infinitesimally  extended  gaps  but 
this  is  not  necessary.  For  gaps  with  finite  width,  the 
same  form  is  valid  with  the  6 -functions  replaced  by 
functions  characterizing  the  gaps. 

We  have  solved  Eq.  (107),  for  the  case  of  two  gaps, 
within  the  framework  of  second  order  perturbation  theory. 
That  is,  we  have  kept  terms  through  ^(a^a^)  .  The  result 
is 


<r  ^  o 


+  f*^e) 
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(108) 

_  X x^i')  ^  i^)  A t»  '  t>-’“)+4i| 

_  ^  Ts - ^'‘-'■xj  V  4>.]}, 

X  ^ 


In  Eq.  (108)  we  have  positioned  the  first  gap  at  x^  =  0 
and  assigned  =  0  . 

This  solution  displays  several  interesting  properties. 
The  solution  in  the  region  0  <  x  <  is  just  that  obtained 
by  expanding  the  large  signal  solution^^^  through  (a^^)  . 

In  the  region  x  >  x^  the  first  two  terms  are  just  the  super¬ 
position  of  the  solutions  which  would  result  if  we  modulated 
the  beam  separately  by  each  gap.  The  last  three  terms  are 
interference  terms  and  occur  because  modulation  exists  on 
the  beam  when  it  enters  the  second  gap.  The  first  two  of 
these  terms  are  what  we  would  expect  if  at  x  =  x^  the  beam 
had  displacement  and  velocity  modulation  respectively  in  the 
amounts  shown  below. 

6.  RADC-TDR-62-199,  IV  Eulerian  Theory,  Section  C, 
Eq.  (IVC-73),  p.  126. 
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The  last  term  is  a  contribution  which  seems  to  mix  the 
arriving  displacement  with  a  second  gap  displacement 
whose  electronic  phase  is  90°  out  of  phase  with  the  ex¬ 
citing  voltage.  The  complexity  of  the  solution  presented 
by  Eq.  (108)  serves  to  demonstrate  the  difficulties  which 
arise  in  a  large  signal  treatment  of  the  multigap  problem. 
It  is  the  interference  terms  which  cause  the  trouble. 
However,  it  is  just  these  terms  which  are  of  great  interest 
in  the  large  signal  situation. 

The  space -charge  field,  velocity  field,  current  and 
charge  density  in  the  modulated  beam  can  be  calculated 
with  the  aid  of  formulas  which  have  been  previously 
derived. 

In  part  2  we  derive  the  integral  equation  (107)  for 
the  multigap  problem. 

In  part  3  we  calculate  the  second  order  perturbation 
theory  solution  given  by  Eq.  (108). 


7.  RADC-TDR-62-199.  IV  Euler ian  Theory,  Section  C, 
part  5,  pp.  128-131. 
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2.  Integral  Equation  for  the  Multigap  Problem 

For  a  beam  modulated  at  N  infinitesimally  extended 
gaps,  the  differential  equation  (8)  for  the  displacement  vector 

(110) 

H 

=  ^  f  X +<r  ~  cd-+ 4>;) 

r»* 

where  x^,  are  the  position  and  phase  associated  with  the 
i^^  gap.  We  define 

M 

F(e,.i  =  I !L‘  ^  c«  4<>t)  ,111) 

Express  (r(x,  6)  in  the  form 

e'’  (112) 

y 

Eq.  ( 1 10)  then  yields 

ill’  =  F,(V.)  (113) 


8.  RADC-TDR-62-199,  IV  Eulerian  Theory,  Section  C, 
Eq.  (IVC-56). 
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(114) 


We  solve  Eq.  (113)  by  Green's  function  methods.  That  is, 


consider  the  function  G, 

X*-)  ^  its  S 


(  F,U')  J/' 


(115) 


(116) 


Thus  G  is  a  solution  to 


(117) 


except  at  the  point  x  =  x'  .  The  solution  for  G  which  satisfies 
the  right  boundary  conditions  is 


re, 

X  'j 

- T - 


(118) 


The  boundary  r'^nditions  have  been  chosen  so  that  the  solution 
(r^(x)  <.<q.  (116)  is  not  affected  by  the  i^  gap  unless 
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X  >  .  With  Eqs.  (118),  (116),  and  ( 1 1 2)  we  obtain 


(119) 


However,  from  Eq.  (114) 


V 


(120) 


and  therefore 

(121) 


y  F,  e- 


Hence,  Eq.  (119)  takes  the  form 


(122) 


Substituting  the  expression  (111)  for  F  in  Eq.  (122)  we  obtain 
the  desired  integral  equation  (107). 
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3.  Second  Order  Perturbation  Theory 
for  the  Two  Gap  Problem _ 


We  expand  <r  in  power  series  in  the  . 


Through  second  order  Eq.  (107)  becomes 


(123) 


Integrating  the  5'  terms  by  parts  and  calculating  the  first 
order  contribution  we  get 

cr  =:  0  )  X  o 


(124) 


<r  ^ 


(125) 


cr  =; 


x>/i 
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We  now  substitute  this  first  order  solution  into  the  right-hand 
side  of  Eq,  (124)  and  perform  the  indicated  integrations.  The 
result  is  given  by  Eq.  (108).  The  factor  1/2  which  multiplies 
the  quantities 


arises  because  the  derivative 

+  ir- 


(127) 


jumps  from  zero  to  (aj^/2)  sin  0  at  x  =  0  and  the  6 -function 
evaluates  this  quantity  at  its  average  value,  (aj/4)  sin  6  . 
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VII.  Large -Signal  Diic  Model  Theory 


1..  Introduction  and  Summary 

We  consider  the  disc  model  from  the  large -signal 
point  of  view.  For  the  infinite  beam  model  the  space- 
charge  forces  are  characterized  by  the  Green's  function 


IC (x-kO  =  i  I)*-)*.'! 


while  for  the  disc  model  we  have 

z 


t  - 


?«■  “  f»’p 

^  5=  be*Kt. 


(128) 


(129) 


We  derive  in  part  2  an  equation  of  motion  for  the 
disc  model  which  is  valid  when  the  signal  level  is  large. 
The  result  is 


-JMx-a'I 

,  0) 


a  -.  •<  Fj  i.A+o',  e'i 


(130) 


-  54  - 


where 


(131) 


is  a  remainder  term  which  vanishes  in  the  case  of  the  infinite 
beam  model.  It  is  conjectured  that  the  remainder  terms  are 
small  and  that  a  good  approximation  for  the  disc  model  may¬ 
be  obtained  by  neglecting  these  terms  in  Eq.  (130). 

In  part  3  we  investigate  the  dispersive  properties  of  the 
waves  characterized  by  the  disc  model  equation  of  motion.  We 
find  that,  while  the  infinite  beam  model  gives  rise  to  two  non- 
dispersive  space -charge  waves,  the  disc  model  on  the  other 
hand  leads  to  four  waves,  all  of  which  are  dispersive.  The 
wave  numbers  in  dimensionless  units  are  as  follows. 


(132) 


Solutions  1,  2  are  space -charge  waves  and  solutions  3,  4  are 
evanescent  waves  which  damp  exponentially. 

The  dispersive  properties  of  the  disc  model  space - 
charge  waves  have  the  effect  of  providing  a  different  space - 
charge  reduction  factor  for  each  harmonic,  r  . 


.  5S  • 


A 


w 


t»‘+  ^T‘ 


> 


(13J) 


This  has  the  effect  of  distorting  the  form  of  the  disturbance 
on  the  beam  as  it  travels  down  the  drift  space. 


2.  Large -Signal  Equation  of  Motion 
for  the  Disc  Model 


(91 

For  one -dimensional  flow'  '  the  equation  of  motion  is 


J  a' 


K  -S  n+e- 


(134) 


In  the  following  we  treat  K(x,  x')  as  a  function  defined  as 
the  limit  of  a  sequence  of  well-behaved  functions.  These 
functions  are  infinitely  differentiable. 


p.  93. 


9.  RADC-TDR-62-199»  IV  Eulerian  Theory,  Section  C, 
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Consider  the  term 


(135) 


The  Taylor  operators  allow  us  to  write  Eq.  (135)  in  the  form 


(136) 


Expanding  integrating  the  n^^  term  n-times 

by  parts  we  get 


(137) 


where 


e.  —  ' 


(138) 


Eq.  (137)  can  be  written 


(139) 


because 
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^  PI  « j  ^^^***:?  j  a'  ■  O 


This  follows  from  the  property 


*  JcU-x'l - o 

t- 


We  expand  the  exponentials  in  Eq.  (139) 


(140) 


(141) 


h{  4x  \A 

'!»»»  » 


*h|  / 


(142) 


The  operation  (d/8x)*^  on  K  is  the  same  as  (*8/8x')^ 
on  K.  Integrating  by  parts,  n>l  times,  Eq.  (142) 
becomes 


I*  > 


(143) 
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Finally  we  get 


•*  k-*-! 


A)«‘ 


Jtl 


(U-*H  +0!  I-! 


(146) 


In  the  case  of  the  infinite  beam  model 


(147) 


(148) 


jL+» 


We  show  that  the  second  term  in  Eq.  (148)  vanishes. 
To  demonstrate  this  we  show  that 

k4l 

7'  C"-')  »-)■“=  ovU='.V- 


(149) 
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We  can  include  the  ternn  m  -  0  elnce  8/8x(l)  s  0 ,  and  hence 
we  wish  to  prove  the  following  lemma. 

k 

^  sff  o  (150) 

h»  »  O  ' 


In  general,  we  define  G(k,  j(  )  by 


S  P 


(151) 


First  we  develop  a  recursion  relation 

k 


icik,/i)  _  7*  (u  .«Or‘r  t.!-- 

■  “4  M  iA*  J 

K*  »  O 


k-i 


•h  s  P^ 


=  k Vy  +  GCl.,i+i) 


Now  consider 


^  L<r-<r)^^o  j  U  — 


(153) 
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and  apply  the  recursion  relation. 

0  -  (154) 

However 

GCW-  1,0)  *  O  k  =  ^3, -  (155) 

so 

GCU.O  =  o  i  ^  =*  -  (156) 

Apply  the  same  reasoning  to  G{k,  1)  .  We  have 

0=  +CCW,0  (1”> 

but 

G  »  O  k  »  -  (158) 

and  hence 

»:  o  J  k  -  -  (159) 

By  induction  we  obtain 

G  Ck,X)  *  O  k  =  i-Hj  i+i, -  (160) 

and  a  special  case  of  this  relation  is  the  desired  lemma, 

GCk,U-0  =  O  (161) 


which  lies,  so  to  speeik,  on  the  edge  of  the  hierarchy  (160). 


.  62  . 


We  have  proved  that  Eq.  (148)  reduce!  to 

I  = 


(162) 


Thus  for  the  infinite  beam  model  Eq.  (134)  becomes 

^  <r  »  -  o(  d)  (163) 


This  equation  was  derived  in  a  more  physical  way  in  a 
previous  report. The  development  above  is  included 
here  because  it  is  needed  for  the  disc -model  discussion 
below. 

For  the  disc -model  we  return  to  Eq.  (146)  which 
applies  to  any  one -dimensional  model.  The  disc -model 
Green's  function^ is 


(164) 


I 

K^i) 


(165) 


10.  RADC -TDR-62-199>  IV  Eulerian  Theory,  Section  C, 

pp.  112-121. 

11.  R.  N.  Wilson,  "Large  Signal  Space -Charge  Theory  of 
Klystron  Bunching,  "  M.  L.  Report  No.  750,  Microwave  Lab., 
Stanford  University,  Calif.,  September  I960,  p.  109. 


-  63  - 


We  substitute  Eq.  (l65)  into  Eq.  (146).  The  result  is 


(166) 


The  last  term  in  Eq.  (166)  can  be  rewritten  as  follows: 

^  (k+O!)''  V(i.k) 


(167) 


Combining  Eq.  (167)  with  Eq.  (166)  the  disc -model  equation 
of  motion  becomes  ^ 

{p/jo 


o<  Ftt  (>‘+^ 

*  ksi 


(168) 


where 


*  lU+«)! 


,169, 
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X 

X 


We  note  that  the  integrand  of  Eq.  (1^9)  vanishes  at 
x'  ,  but  the  peak  of  the  weighting  function,  W,  lies  at 
x'-. 


(170) 


4*0 


1 


(171) 


The  quantity 


is  a  weighted  average  of  the  quantity 


U  +  O! 

about  the  value  x'  =  x  ,  with  weighting  function  W  .  (Rh 
vanishes  in  the  two  limits 


(172) 


V,  ^  o  ,  yg  — ► 


] 


(173) 


where  b  is  the  beam  diameter.  We  conjecture  that 


y  ^  rf?i, 

k*i 


(174) 


is  small  compared  to  the  other  terms  in  Eq.  (I6d).  In  this 
case  a  good  approximation  to  the  disc -model  equation  of 
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motion  is 

(175) 


.f_i£ 


s=  -  0^  ©■) 


3.  Dispersion  Characteristics  of  the  Approximate 
Large -Signal  Disc  Model _ 

To  determine  the  type  of  waves  which  the  disc- 
model  supports,  we  deal  with  the  homogeneous  version 
of  Eq.  (175)  and  express  (r(x.  0)  in  the  form 

<rC^e)*  (rlo,  p) 

Also  we  express  the  kernel  function  in  the  integral  term  of 
Eq.  (175)  in  terms  of  its  Fourier  transform. 


-I-® 


(177) 
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Upon  substitution  into  the  homogeneous  form  of  Eq.  (175)  the 
dispersion  relation  which  results  is 


(Jeb 


(178) 


If  we  solve  for  P(u)  we  obtain  four  roots.  Through  leading 
order  in  these  are 


1'^ 


(uUrJ  J 


(179) 


The  waves  |3,  ,  are  space -charge  waves  and  |3_  .  are 
evanescent  waves. Thus  whereas  there  were  only  two 
non -dispersive  space -charge  waves  for  the  infinite  beeun 
model,  we  now  have,  in  addition,  two  evanescent  waves. 
Also  all  four  of  the  waves  characterized  by  Eq.  (179)  are 
dispersive  and  this  makes  the  situation  more  complicated. 
For  modulation  by  an  infinitesimally  extended  gap  a 

(13) 

first  order  solution  has  previously  been  obtained.  We 
shall  here  neglect  the  evanescent  waves  and  match  boundary 


12.  R.  N.  Wilson,  M.L.  Report  No.  750,  pp.  112-116. 

13.  I]^,  pp.  112-113. 
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conditions  in  the  same  way  as  for  the  infinite  beam 
model.  We  let 


(Tio^e)  *=  o 


t-o 


(180) 


and  approximate  the  second  of  Eqs.  (180)  by 


a 


(181) 


We  obtain  the  following  solution. 


0  j  <  o 


er  (n  fl)  = 


=^1 

r  S*  ~\ 

>***(|H1V*'* 

L  -J 

(d’X)^  K  >  O 


This  solution  agrees  with  the  solution  obtained  before, 
at  points  down  the  drift  space  such  that  the  evanescent 
waves  (exponentially  damped)  are  negligible. 

This  suggests  that  a  good  approximate  solution  is 
provided  by  neglecting  the  evanescent  waves  when  we  fit 


(182) 


14.  RADC -TDR -62-199,  IV  Euler ian  Theory,  Section  C, 
p.  124. 
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the  boundary  conditions.  The  resulting  solution  is  meaning¬ 
ful  only  at  distances  from  the  gap  which  are  large  enough  so 
the  evanescent  waves  can  be  neglected. 

We  approximate  Eq.  (180)  by 

<r(0,  e)  =  O 

ie  i  ^  ^ 

'  *=  4  A^6  - 

The  d-c  part  of  Eq.  (183)  gives  rise  only  to  evanescent  waves. 
The  solution  is 


(T  (X|0)  “ 


0  ^  X  <  O 

{%) 

0  >  0 


TTTTh''*  j 


tx 

L  fvt/vN 


(184) 


Eq.  (184)  shows  that  there  is  a  different  space-charge  reduction 
factor  for  each  harmonic. 


^  = 


U)  yr 
U) 


Vy 


u> 


Vt‘+ 


y 


(185) 
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In  the  case  of  the  infinite  beam,  =  1  ,  we  found  that  the 
disturbance  on  the  beam,  caused  by  the  modulation,  traveled 
down  the  drift  space  without  change  in  form  but  changed  in 
amplitude  by  the  multiplicative  factor 

- - -  (186) 

This  is  no  longer  true  with  the  disc  model  as  shown  by  Eq. 

(184).  Each  harmonic  has  a  different  space -charge  multi¬ 
plicative  factor 


4  •  ^ 

H  A 


(187) 


and  the  reason  for  this  behavior  lies  in  the  dispersive 
properties  of  the  space -charge  waves. 


-  70  - 


VIII.  Statistical  Machaniol  Formulation 
of  the  Infinite  Beam  Model 


1.  Introduction  and  Summary 

We  formulate  the  infinite  beam  problem  within  the  frame¬ 
work  of  statistical  mechanics.  We  start  with  the  statement  of 
Liouville's  theorem  for  the  many  particle  density  function  and 
derive  from  this  an  equation  of  motion  for  the  electron  beam. 
The  result  is 


(189) 


=  d-c  beam  velocity 
e(n  +  6n)  =  electron  charge  density 


(190) 


The  term  on  the  right  is  a  diffusion  term  which  appears  to 
beconne  large  at  and  near  crossover. 
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We  conclude  that  the  infinite  beam  model,  which 
pictures  the  beam  as  a  set  of  charge  planes,  ceases  to 
be  valid  near  crossover. 


2.  Derivation  of  the  Hydrodynamical 
Equations  of  Motion _ 

We  start  with  the  many  particle  distribution 
function,  ^ 


(191) 


for  two  species  of  particles,  electrons  (q^,  p^)  and  ions 
(Q^,  P^)  .  According  to  Liouville's  theorem  the  total  time 
derivative  of  W  vanishes. 


The  positions  q^,  Q.  and  momenta  p^,  P.  satisfy 

Ti” 


15.  R.  C.  Tolman,  The  Principles  of  Statistical  Mechanics, 
Oxford  Univ.  Press,  New  York,  N.Y. ,  1st  e^tion,  pp.  48-52. 
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We  define  the  single  particle  electron  distribution  function 


by 


"  j'W  Cl;, a;,  f..,?.) 


and  noi'malize  the  density  function  W  so  that 


(195) 


(196) 


Eq.  (195)  gives  the  electron  particle  density  and  Eq.  (196) 
gives  the  total  number  of  electrons. 

Integrate  Eq.  (192)  over  all  coordinates  except  q^,  p^  . 
The  result  is 


n  -f.  :  O 


4  1 X; 


4  -• 


^  ifi 

i  =  • 


(197) 


All  other  terms  vanish  since,  for  example 


(198) 
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However,  the  last  term  vanishes  because  W  is  assumed 
to  vanish  at  the  boundary  of  integration  in  the  phase  space. 
We  substitute  for  p^  in  Eq.  (197),  from  Eq.  (193)  and 

obtain 


O 


(199) 


as  the  differential  equation  satisfied  by  the  single  particle 
distribution  function. 

Note  that  we  have  not  treated  the  electromagnetic  field 
statistically.  This  implies  an  assumption  as  to  the  macro- 
scopic  reproducibility  of  E,  H  . 

For  our  beam  model  we  write  fa>(x,  p)  in  the  form 


(  200) 


where 


“o  =  -  P, 


p  =  mv 
*^o  o 


(201) 


Pq  b«3ini  electron  momentum  and  we  have  assumed 

that  the  d-c  beam  is  uniform  with  electron  charge  density  en . 
In  what  follows  wo  shall  neglect  magnetic  effects.  Eq.  (199) 


becomes 


-  74  - 


u 


+ 


^  irj  + 


^  Pm. 


O 

(202) 


To  obtain  hydrodynamical  equations  of  motion  we  define 


(203) 


th 

and  take  the  zero  and  first  moments  of  Eq.  (202)  with  respect 
to  p  .  We  make  use  of  the  following  definitions  for  average 
values . 

p  =  e  CuA’p  =  ^  j  (u.  +  Mljp 

-T  e v\  +  e 


(204) 


(205) 


(206) 


(207) 
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A’p=  Si*'?.  +  {iu  p«l’p} 


r  4,  +  a 


(208) 


First  we  multiply  Eq,  (202)  by  e  and  integrate  over 
th 

p  .  The  resulting  zero  moment  equation  is 


>1*^  "  ° 


(209) 


since  the  other  two  integrals  vanish.  Eq.  (209)  is  the  equation 
of  continuity  of  charge  for  the  electrons. 

To  calculate  the  first  moment  equation  we  multiply 
Eq.  (202)  by  ep^/m  and  integrate  over  p. 


"Hu..  + 


(210) 


With  Eq.  (203)  the  integral  term  in  Eq.  (210)  can  be  rewritten 


^  p«>i.  f*v  +  Q  pMc  A  Pi>  +  fi/6  A  p>.  ) 

..j.,  +  i*.)  A  pM  A  Pjj  ^  ^ 


(211) 


We  rewrite  the  first  term  in  Eq.  (210)  as  follows.  From 
Eq.  (208) 

i.a  =  "  ft  fvA’p  -  ® 


(212) 


(K+Jh.)Apt> 
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Hence 

u 


it 


(213) 


but  using  Eqs.  (209),  (207)  aud  (212)  we  can  write 


u 


(214) 


Combining  Eq.  (213)  and  Eq,  (214)  we  get 

^  =  -  -^  ('po.  P*- 


(215) 


Substituting  Eqs.  (211)  and  (215)  into  Eq.  (210)  and  re¬ 
arranging  terms  we  finally  obtain 

+  “i:  ^Pv)  -  e 


(216) 


^  (k  4  Tk.)  L 

which  is  the  hydrodynamical  equation  of  motion.  The  left- 
hand  side  is  the  well-known  equation  of  motion  for  the  beam 
in  Eulerian  coordinates 
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-  e  e  =  0 


(217) 


r  refers  to  the  electron  position  in  the  presence  of  the 
modulation.  ^ 

The  term  on  the  right-hand  side  of  Eq.  (216)  is  a  heat 
diffusion  term. 


3.  Remarks  on  the  Large -Signal 
Infinite  Beam  Model _ 

We  specialize  Eq.  (2l6)  to  the  case  where  mAv  = 

Ap  =  1  Ap(z,t)  and  E  =  1  E(z,t).  The  result  is 
z  z 

p/ii  + 

_ ^  (218) 

We  have  solved  Eq.  (218),  neglecting  the  right-hand  side, 
within  the  framework  of  the  infinite  beam  model.  The 
solution  was  valid  up  to  the  point  where  crossover  sets  in. 

Eq.  (218),  on  the  other  hand,  is  a  more  general  equation. 

It  was  pointed  out  in  our  discussion  of  the  infinite  beam 
model  that  at  crossover  the  charge  density  becomes  infinite. 


pp. 


16.  RADC -TDR -62-199,  IV  Eulerian  Theory,  Section  A, 
62-63. 

17.  Ibid,  Section  C,  pp.  91-139. 
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That  is,  in  the  quasi -Eulerian  formalism  of  our  previous  work 
the  charge  density,  at  the  actual  position  of  the  electron  is 
given  by 

'  J>0 

and  1  +  dcr/dx^  passes  through  zero  at  crossover. 

It  was  mentioned  in  connection  with  our  previous  analysis 
that  one  might  question  the  validity  of  the  infinite  beam  model 
at  and  beyond  crossover.  The  infinite  beam  model  assumes 
that  one  can  describe  an  infinitely  extended  ion  neutralized 
electron  plasma  by  a  set  of  charge  sheets  constrained  to  move 
in  one  (z)  direction.  The  derivation  of  Eq.  (218)  has  avoided 
the  charge  sheet  picture  and  enables  us  to  make  statements 
regarding  the  validity  of  the  infinite  beam  model.  Eq,  (218) 
states  that  the  infinite  beam  model  is  a  good  approximation 
when  the  right-hand  term  is  negligible.  However,  let  us  look 
at  this  term  near  crossover.  We  rewrite  it  in  the  form 


where 


?  =  - 


(221) 
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At  crossover,  according  to  infinite  beam  model  theory, 
n  +  5n  passes  to  infinity  and  it  would  take  a  peculiar  can¬ 
cellation  to  keep  the  expression  (220)  from  becoming  very 
large.  On  the  other  hand,  the  terms  on  the  left  side  of 
Eq.  (218)  add  up  to  zero  within  the  framework  of  the  infinite 
beam  model.  We  therefore  conjecture  that  the  infinite  beam 
model  ceases  to  be  valid  at  and  near  crossover. 
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IX,  Infinite  Beam  Model  With 
Phenomenological  Losses 

1 .  Introduction  and  Summary 

In  our  statistical  mechanical  treatment  of  the  infinite 
beam  model  we  found  that  the  hydrodynamical  equation  of  motion 
is  an  approximate  equation  which  neglects  heat  diffusion  type 
terms.  We  attempt  here  to  include  such  effects  in  an  approx¬ 
imate  way  by  the  introduction  of  a  collision  loss  term  into  the 
equation  of  motion.  This  device  is  common  in  plasma  theory 
and  is  used  to  represent  the  effect  of  collisions  of  the  electrons 
with  neutrals.  In  the  present  situation  this  term  should  pro¬ 
vide  a  step  in  the  direction  of  a  more  realistic  beam  model. 

In  part  2  we  derive  an  equation  of  motion  in  the  presence 
of  losses.  The  result  is 

(222) 

The  factor  2^11  represents  the  number  of  "collisions"  per 
r-f  cycle  of  the  modulating  voltage. 

In  part  3  we  are  able  to  obtain  a  large -signal  solution 
to  Eq.  (222)  for  the  case  of  modulation  at  an  infinitesimally 
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extended  buncher  gap.  This  solution  appears  in  closed  form 
and  contains  complete  information  about  harmonics  of  all 
orders.  The  solution  is 


<r(* 


,e,, 

!>>  PfJ 


”'3 

(223) 


where 


(224) 


In  part  4  we  remark  on  crossover  criteria  in  the 
presence  of  the  damping  term.  For  no  losses,  crossover 
occurs  when  the  depth  of  modulation,  6,  is  approximately 
equal  to  2.0.  Employing  the  same  criterion  for  the  lossy 
case  we  find  that  crossover  is  inhibited  by  the  damping. 
The  critical  depth  of  modulation  moves  from  2.  0  in  the 
case  of  no  damping  to  the  value  5.43  in  the  case  of  critical 
damping. 
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2.  Equation  of  Motion  With  Phenomenological 
Collision  Losses 


The  infinite  beam  model  hydrodynamical  equation  of 
motion  in  the  absence  of  losses  is 


where 


X  = 

ieL  ^ 

(r  = 

O  •a 

0  = 

X  = 

(226) 


(227) 


z  is  distance  down  the  drift  space,  and  ^  is  the  displacement 
caused  by  the  presence  of  the  modulation. 

In  our  statistical  mechanical  treatment  of  the  infinite 
beam  model  we  found  that  Eq.  (226)  is  an  approximate  equation 
which  neglects  heat  diffusion  terms.  We  shall  here  attempt  to 
make  a  correction  for  these  extra  terms  by  introducing  into 
Eq.  (226)  a  phenomenological  collision  loss  term.  We  con¬ 
sider  that  the  motion  is  degraded  by  collision  losses  charac¬ 
terized  by  a  collision  frequency  v  which  is  taken  independent 
of  velocity.  Thus  we  enter  a  term  in  the  equation  of  motion  of 
the  form 


w\  f  ^  ■$  Uii) 


(228) 


Furthermore  we  assume  that  the  d-c  beam  is  produced  in 
a  way  that  collision  losses  in  the  unmodulated  beam  are 
accounted  for.  so  that  we  are  left  with  a  beam  traveling  at 
constant  velocity  v^  along  the  z-axis. 

In  dimensionless  units  Eq.  (226)  now  becomes 

^  p/i<9  +  ^^1)*  ]  -f  n 

where 

i’!’S=''A  J 

ZtTt]  is  the  number  of  collisions  per  r-f  cycle  of  the  gap 
modulating  field. 

The  device  of  introducing  such  a  phenomenological 
collision  damping  term  is  familiar  in  plasma  theory  to 
account  for  collisions  of  electrons  with  neutral  atoms  and 
molecules.  In  our  present  application  this  term  probably 
provides  a  step  in  the  direction  of  a  more  realistic  beam 
model. 

3.  Large -Signal  Solution  of  the  Equation  of  Motion 

We  consider  the  problem  of  modulation  by  an  in¬ 
finitesimally  extended  gap  located  at  x  =  0  .  The  dis¬ 
placement  ff  is  expressed  in  the  form 
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(229) 


(230) 
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^(>‘.<9)  =  ^7v(k)  e 

Y 

(231) 

Substituting  this  expression  into  the  homogeneous  form  of 

Eq.  (229)  we  obtain  the  following  equation  for  the  «r^(x)  , 

- 

(232) 

We  express  <r^  by  a  Fourier  integral 

TV  =  ^  TV((J)  e 

(233) 

This  results  in  the  following  dispersion  relation. 

• 

(>*+  ipn  -A‘  =  o  'A 

(234) 

5  iX 

Hence  (r  takes  the  form 

) 

4<rv» 

(235) 
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We  treat  the  gap  as  a  boundary  condition  on  r . 

-  o  I 

“IT  “7^ j  = 

)t  S  O  ^ 

The  first  of  Eqs.  (236)  yields 

-  cr,,  (237) 


or 


(^0)  >ua.[;cA).;jp7C']  ^  ^ 


(238) 


Making  use  of  the  second  condition  (236)  we  finally  obtain 


AJTF 


(239) 


(18) 

This  solution  reduces  properly  to  the  lossless  case'  '  when 


=  -2!L  - ^  o 

iX 


(240) 


18.  RADC -TDR -62-199,  IV  Eulerian  Theory,  Section  C, 

p.  126. 
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Eq.  (239)  exhibits  a  new  feature,  namely,  that  the  displacement, 
(T  ,  caused  by  the  modulation  decays  exponentially  with  x  .  The 
1/e  point  is  given  by 


= 


rx 


(241) 


4.  Crossover  Criteria 


It  is  interesting  to  investigate  the  effect  of  the  loss  term, 
upon  crossover  criteria.  Crossover^ occurs  when 

)+  O  (242) 


From  Eq.  (239) 
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(243) 


19.  RADC -TDR-62 -199,  IV  Eulerian  Theory,  Section  C, 
pp.  137-139. 
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To  obtain  some  feeling  for  the  crossover  conditions  we 
consider  an  electron  such  that 


6-y  ’ 


(244) 


For  this  phase 


^-fCAy) 


(245) 


Crossover  does  not  occur  here  if 


I  (246) 


The  maximum  of  I(x)  occurs  at 


! —  W 


(247) 


Hence 


XjT^  f---,  =  yTi^ 


I  (  ^  j 
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(249) 


We  define  the  depth  of  modulation,  6,  by 


(250) 
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The  condition  (249)  becomes 

4-  < 


(251) 


where 


c(g)  = 


(252) 


C(^)  is  plotted  below.  4  =  0  corresponds  to  no  damping  and 
the  crossover  criterion  (251)  is 


(253) 


4=1  corresponds  to  critical  damping  and  in  this  case  the 
crossover  criterion  becomes 


S<  2«e- 


(254) 
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We  see,  therefore,  that  the  presence  of  the  collision 
term  tends  to  inhibit  crossover:  i.e. ,  the  depth  of  modulation 


(255) 


is  allowed  to  have  larger  values  in  the  presence  of  damping 
before  crossover  sets  in. 
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X.  Remarks  on  the  Possible  Directions 
of  Future  Work 


We  remark  briefly  on  the  directions  which  future  work 
in  beam  theory  might  take. 

1.  The  vector  normal  mode  formalism ,  although 
complicated,  contains  a  fairly  complete  picture  of  the  inter¬ 
action  of  an  electron  beam  with  some  attendant  electromagnetic 
structure.  It  provides  an  accurate  stai'ting  point  for  numerical 
analysis  of  specific  devices  in  the  large  signal  regime.  The 
analysis  could  profitably  be  extended  to  the  case  of  magnetically 
focused  beams  which  allow  radial  motion,  in  particular  to  the 
Brillouin  beam. 

2.  We  have  made  sizable  progress  in  the  large -signal 
description  of  a  modulated  electron  beam  within  the  framework 
of  the  infinite  beam  model.  An  accurate  description  has  been 
obtained  for  the  case  of  one  infinitesimally  extended  modulating 
gap  and  the  analysis  avoids  the  use  of  perturbation  theory.  Some 
work  has  been  accomplished  on  the  two -gap  problem  with  the  use 
of  second  order  perturbation  theory.  The  multigap  problem 
should  be  analyzed  in  the  large -signal  regime  by  methods  which 
avoid  perturbation  theory. 

20.  RADC -TDR -62-199,  IV  Eulerian  Theory,  Section  B 
pp.  71-91. 


Some  work  has  been  done  which  extends  the  large - 
signal  infinite  beam  analysis  to  the  case  of  finite  geometry, 
specifically,  the  disc -model.  This  is  a  fertile  area  for 
future  endeavor. 

3,  Our  recent  work  with  the  infinite  beam  model 
indicates  the  possible  breakdown  of  the  hydrodynamic 
theory  at  and  near  crossover.  Thus  if  one  wishes  to  push 
large -signal  beam  theory  into  the  region  where  crossover 
is  imminent,  the  problem  should  be  formulated  from  the 
point  of  view  of  statistical  mechanics.  The  reason  for  this 
is  that  terms  which  are  neglected  in  the  hydrodynamic 
approximation  appear  to  be  significant  in  the  region  near 
crossover. 
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